Precision sparticle spectroscopy in the inclusive same-sign dilepton channel at LHC 



O 

D 



Oh' 
I 

Oh' 

CD 



> 
O 
O 
m 

G\ 
O 
0^ 
O 



X 



Konstantin T. Matchev^, Filip Moortgat^, Luc Pape^, and Myeonghun Park^ 
^Physics Department, University of Florida, Gainesville, FL 32611, USA and 
'^ETHZ, Zurich, Switzerland 
(Dated: 27 July, 2010) 

The inclusive same-sign dilepton channel is already recognized as a promising discovery signature 
for supersymmetry in the early days of the LHC. We point out that it can also be used for precision 
measurements of sparticle masses after the initial discovery. As an illustration, we consider the 
LM6 CMS study point in minimal supergravity, where the same-sign leptons most often result from 
chargino decays to sneutrinos. We discuss three different techniques for determining the chargino 
and sneutrino masses in an inclusive manner, i.e. using only the two well measured lepton momenta, 
while treating all other upstream objects in the event as a single entity of total transverse momentum 
Pt- This approach takes full advantage of the large production rates of colored superpartners, but 
does not rely on the poorly measured hadronic jets, and avoids any jet combinatorics problems. We 
discuss the anticipated precision of our methods in the early LHC data. 

PACS numbers: 14.80.Ly,12.60. Jv,11.80.Cr 



A long standing problem in hadron collider phe- 
nomenology has been the determination of the absolute 
mass scale of new particles in events with missing en- 
ergy. The prototypical example of this sort is provided 
by any model of low-energy supersymmetry (SUSY) with 
conserved i?-parity, in which the lightest superpartner 
(LSP), typically the lightest neutralino Xii is a neutral, 
weakly interacting particle of a priori unknown mass [l| . 
Astrophysics also adds credence to such scenarios, since 
the LSP is a potential dark matter candidate, whose 
relic abundance is typically in the right ballpark R- 
parity conservation guarantees that every event contains 
(at least) two invisible particles, whose energies and mo- 
menta are not measured, making the full reconstruction 
of such events a very challenging task. 

Recently, several solutions to this problem at the Large 
Hadron Collider (LHC) have been proposed. Most of 
them rely on exclusive channels Q, where a sufficiently 
long decay chain can be properly identified. Unfortu- 
nately, this almost inevitably requires the use of hadronic 
jets in some form in the analysis - in most SUSY mod- 
els, the main LHC signal is due to the strong production 
of colored superpartners, whose cascade decays to the 
neutral LSP necessarily involve hadronic jets. For many 
reasons, jets are notoriously difficult to deal with, espe- 
cially in a hadron collider environment. Because of the 
high jet multiplicity in SUSY signal events, any jet-based 
analysis is bound to face a severe combinatorial problem 
and is unlikely to achieve any good precision. Thus it is 
imperative to have alternative methods which avoid the 
direct use of jets and instead rely only on the well mea- 
sured momenta of any (isolated) leptons in the event. 

In this letter, we describe three such methods, which 
are free of the jet combinatorial problem. For illustration, 
we shall use the standard example of i?-parity conserv- 
ing supersymmetry with a Xi LSP. Its collider signatures 
have been extensively studied, and typically involve jets, 
leptons and missing transverse energy [J. Among those, 
the inclusive same-sign dilepton channel has already been 



TABLE I: Selected sparticle masses (in GeV) at point LM6. 
We list the average ql mass Mgj^ = ^{Mn,^ + M^^). 
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identified as a unique opportunity for an early SUSY dis- 
covery at the LHC 0, Q . The two leptons of the same 
charge can be easily triggered on, and provide a good 
handle for suppressing the SM background. In our anal- 
ysis we use the LM6 CMS study point 4] , whose relevant 
mass spectrum is given in Table HI At point LM6, sig- 
nal events with two isolated same-sign leptons typically 
arise from the SUSY event topology in Fig. [T] Consider 
the inclusive production of same-sign charginos, which 
decay leptonically as shown in the yellow-shaded box in 
the figure. The resulting sneutrino could be the LSP 
itself, or, as in the case of LM6, may further decay invis- 
ibly to a neutrino u and the true LSP Xi- Such same- 
sign chargino pairs typically result from squark decays, 
as indicated in Fig. [TJ In turn, the squarks may be pro- 
duced directly through a t-channel giuino exchange, or 
indirectly in giuino decays. Note that the two same-sign 
leptons in Fig. [1] are accompanied by a number of up- 
stream objects (typically jets) which may originate from 
various sources, e.g. initial state radiation, squark de- 
cays, or decays of even heavier particles up the decay 
chain. In order to stay clear of jet combinatorial issues, 
we shall adopt a fully inclusive approach to the same-sign 
dilepton signature, by treating all the upstream objects 
within the black rectangular frame in Fig. [T] as a single 
entity of total transverse momentum Pt- 

Given this very general setup, we now pose the follow- 
ing question: assuming that a SUSY discovery is made 
in the inclusive same-sign dilepton channel, is it possible 
to measure the individual sparticle masses Mp and Mc 
involved in the leptonic decays of Fig. [1] using only the 
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FIG. 1: The typical SUSY event topology producing two iso- 
lated same-sign leptons at point LM6 (see text for details). 
The diagram for a pair of negatively charged leptons l~ l~ is 
analogous. 

transverse momenta of the two leptons p^j} and p^^^ , and 
the total upstream transverse momentum P^? Although 
it may appear that those three vectors do not provide a 
lot of information to go on, we shall show that this is 
possible. We discuss three different approaches. 

Method I. Let us concentrate directly on the observed 
lepton momenta p^j, . Consider the two collinear momen- 
tum configurations illustrated in Fig. [2] and defined as fol- 
lows. In each configuration, the lepton momenta are the 

same: p^j! — Pg^; and then they can be either parallel 
or anti-parallel to the measured upstream Pt'. 



Pit = PeT tt Pt 
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to refer to the parallel (anti-parallel) configuration: s = 
^^\Pt) = cos{p^^,Pt)- Now let us measure the 
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maximum lepton momentum in each configuration: 
Pit{sPt)= max _ _ Ipir}- 
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Observe that both piT^+PT) and piT{—PT) can be di- 
rectly measured from the lepton pT distributions. For 
example, construct a 2D scatter plot {x,y} of 



X = cos{pIj} +pI^' ,Pt), 
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with the cut jp^j. — p^j. | < e ('^ 0), and take the limit 

Pit{sPt) = lim . (5) 

Armed with the two measurements p^t(+^t) and 
Pit{—Pt), we can now directly solve for the masses Mp 
and Ale- The formula for ptT{sPT) is 
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(^4M2 + (sPr)2 - sPt) ■ (6) 



Inverting ([6]), we get 
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FIG. 2: The two special momentum configurations defined in 
eqs. (HEl). 



thus fixing the absolute mass scale in the problem. Once 
the parent mass Mp is known, the child mass Mc is 



Mr 




Pit{-Pt) - Pit^+Pt) 



(8) 



Thus we found the true sparticle masses Mp and Mc 
directly in terms of the measured lepton momenta 
Pit{'^Pt) and upstream momentum Pt- Note that the 
choice of the value for Pt in eqs. ^ and ([8]) is arbitrary, 
which can be used to our advantage, e.g. to select the 
most populated Pt bin, minimizing the statistical error. 

Method II. In our previous method, the lepton mo- 
menta pit{^Pt) were measured directly from the data 
as implied by eq. Alternatively, we can obtain them 
indirectly from the endpoint of the Cambridge Mt2 vari- 
able 0. To be more precise, we apply the "subsystem" 
Mt2 variable introduced in 7] to the purely leptonic sub- 
system in the yellow-shaded box of Fig. [T] Following the 
generic notation of Rcf. 7], we denote the input (test) 
mass of the sneutrino child as Mc- The subsystem Mt2 
variable is now defined as follows. First form the trans- 
verse mass Mt for each (chargino) parent 



M2 + 2( |p«|VM2-H|p; 



^cT 



' Pit ■ PcT 



in terms of the assumed test mass Mc and transverse 
momentum p^,^ for each (sneutrino) child. Just like the 
traditional Mt2 the leptonic subsystem Mt2 vari- 
able is defined through a minimization procedure over 
all possible partitions of the unknown children momenta 
p^J^ , consistent with transverse momentum conservation 

Mt2{Mc, Pt.pIt) = min {max |il/4^\ M^^^ 1 1 . (9) 
The Mx2 distribution has an upper kinematic endpoint 
M:^^''{Mc,Pt)= max {Mt2{Mc,Pt,P^^)] . (10) 

all events L J 

which can be experimentally measured and subsequently 
interpreted as the corresponding parent mass Mp 



Mp{Mc, Pt) = M^riMc, Pt) , 



(11) 



Mr. 



\/peT{-PT)PiT{+PT) 

Pit{~Pt) -Pit{+Pt) 



PT: 



(7) 



providing one functional relationship among Mp and Mc, 
but leaving the individual masses still to be determined. 
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For us the importance of the Mt2 variable ^ is that 
the monientuni configurations in Fig. [5] are precisely the 
ones which determine its endpoint M^2^ . The complete 
analytical dependence of the Mt2 endpoint Mp{Mc, Pt) 
on both of its arguments Mc and Pt is now known Q : 



MpiM^Pr) = 
where 

Mp{M,,sPt) 



MpiMc+Pr), if Mc<M„ 
Mp{M„-Pt), if Mc>M„ 
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Thus we can alternatively obtain the sparticle masses by 
measuring just two Mt2 kinematic endpoints, with arbi- 
trary choices for the test mass Mc and the upstream Pt- 
For concreteness, let us pick some fixed and P^, form 
the corresponding Mt2 distribution © and measure its 
endpoint M^, also making a note of the configuration s': 



r A>' D' 1 measure f ,\ 

Now perform a second such measurement 

f tT/T" d"1 -measure ( ,->// i,\ 



(14) 



(15) 



By inverting (|T3|). these two measurements allow the ex- 
perimental determination of 
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and similarly for pit{s" P!P). Now taking the ratio 
_ piT{s'P!r) _ 



4M2 + (s'P^)2 - s'P^ 
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where in the second step we used eq. ([6]), we can solve 
()17ll for the true parent mass Mp in terms of measured 
quantities: 



Mp = 



-rs'P^s"P/^ 



(1 



s"P^' 
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(18) 



and then find the true child mass A/c from (El) as 



M, = M„ 



with Mp already given by (IT51) . Note than in this method, 
the values of M^, Af", P^ and P^' can be chosen at will, 
allowing for repeated measurements of Alp and Mc- 
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FIG. 3: M^2 ^ versus the test mass Mc, as obtained in our 
simulations (data points) from a sample with Pt = 420 ± 50 
GeV, or theoretically from eq. (I12p (blue solid line), as well 
as their difference (lower panel). 



Method III. The third and final method for extract- 
ing the two masses Mp and Mc will make use of the 
celebrated "kink" in the Mt2 endpoint function p2)) 
01. Since Afp(Mc,+PT) and MpiMc-Pr) have differ- 
ent slopes at the crossover point Mc — Mc, the function 
Mp{Mc,Pt) has a slope discontinuity precisely at the 
correct value Mc of the child mass, providing an alter- 
native measurement of the absolute mass scale The 
procedure is illustrated in Fig. |3] for the LM6 study point 
of Table HI The blue solid line shows the theoretically 
expected shape from eq. ([T2|) . for Pr = 420 GeV, which 
is roughly the mean of the Pt distribution at point LM6. 
In the LM6 case the kink is very mild, only 3.3° 0. 

In order to test the precision of the three methods, we 
perform event simulations using the PYTHIA event gen- 
erator |9|] and PGS detector simulation ^Ifl]. We consider 
the LHC at its nominal energy of 14TeV and 100 fb"^ of 
data. To ensure discovery, we use standard CMS cuts as 
follows [1, HH: exactly two isolated leptons withp^ > 10 
GeV, at least three jets with pt > (175, 130, 55) GeV, 
^T > 200 GeV and a veto on tau jets. With those cuts, 
in the dimuon channel alone, the remaining SM back- 
ground cross-section is rather negligible (0.15 fb), while 
the SUSY signal is 14 fb, already leading to a 22a discov- 
ery with just 10 fb^^ of data [1, [HI. In order to compare 
to the theoretical result in Fig. |3l we select a ±50 GeV 
Pt bin around Pt = 420 GeV and construct a series of 
Mt2 distributions, for different input values of Mc- For 
each case, we include all SM and SUSY combinatorial 
backgrounds, and extract the M^2^ endpoint by a lin- 
ear unbinned maximum likelihood fit, obtaining the data 
points shown in Fig. |31 We see that the Mt2 endpoint 
can be determined rather well {5Mp ^ 3 GeV), but only 
on the right branch Mc > Mc- In contrast, the Mt2 end- 
points on the left branch Mc < Mc are considerably un- 
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FIG. 4: Scaling factors relating the error SMp in the extrac- 
tion of the Mt2 endpoint to the resulting uncertainties 5Mf 
and 5Mc on the parent and child masses calculated from (|18[ 
and (|19p . as a function of the true input masses Mc and Mp 



derestimated, washing out the expected kink. There are 
two separate reasons behind this effect. Recall that the 
Mt2 endpoint on the left branch is obtained in the config- 
uration s = +1 of Fig. [21 which requires the lepton to be 
emitted in the backward direction. As a result, the par- 
ent boost favors configurations with s ~ — 1 over s ~ 
Another consequence is that leptons with s ~ -f 1 are 
softer and more easily rejected by the offline pr cuts. We 
conclude that M^2^ measurements on the left branch are 
in general not very reliable, and tend to jeopardize the 
traditional kink method. For example, using Method III 
to fit the data in Fig. [3] (green dotted line), we find best 
fit values of only Mp(fit) — 212 GeV and Mc^fu) — 188 
GeV. Method I has a similar problem, since pit{+Pt) is 
measured from events in the s = +1 configuration. Using 
the Mp measurements from Fig. [3] at Mc = and Mc = 1 
TeV, we find from eq. ^ that ^^(+420 GeV) = 8.8 
GeV and ^^(-420 GeV) = 50.6 GeV (compare to the 



nominal values of 14.8 GeV and 53.6 GeV, correspond- 
ingly). The resulting mass determination via eqs. (j7l8p 
is Mp(^fit) = 212 GeV and M^^^n-^ = 190 GeV. We see 
that in both Method I and Method III, the masses are 
underestimated due to the systematic underestimation of 
the left M^2^ branch in Fig. [3l It is therefore of great 
interest to have an alternative method, which relies on 
the right M^2 -^ branch alone. 

This is where the available freedom in Method II comes 
into play, since both test masses M^ and M" can be cho- 
sen on the right branch. Taking = 350 ± 50 GeV and 
P^' — 500 ± 50 GeV and repeating our earlier analysis, 
we find that 5Mp on the right branch is still on the order 
of 3 GeV, as in Fig. [3l The resulting error 6Mp (SMc) 
on the measured parent (child) mass can be easily prop- 
agated from eqs. (|18ll9p . The two ratios SMp/ SMp and 
SMc/ SMp are shown in Fig. 21 where for concreteness we 
have taken M^ = M'J = 1000 GeV. Fig. g] reveals that 
the LM6 input values of Mc and Mp are rather unlucky, 
since the error SMp on the Mt2 endpoint is then ampli- 
fied by a factor of almost 70. However, if Mc and Mp hap- 
pened to be different, with the rest of the spectrum the 
same, the precision quickly improves. For example, with 
SMp — ±3 GeV, the masses can be determined to within 
±30 GeV (±75 GeV) within the yellow (orange) region. 
One should keep in mind that the dominant uncertainty 
on SMp is due to the SUSY combinatorial background. 
We have verified that in the absence of such combinato- 
rial background, SMp <, 1 GeV and the typical precision 
on Mp and Mc from Fig. [4] is then at the level of 10%. 

In conclusion, we considered the inclusive same-sign 
dilepton channel in SUSY, which so far has only been 
used for discovery, but not for mass measurements. We 
demonstrated that it allows a separate determination of 
the chargino and sneutrino masses. We discussed three 
different methods, which rely exclusively on the well mea- 
sured lepton momenta. The methods are completely gen- 
eral and inclusive, and can be applied to other SUSY 
topologies and to non-SUSY scenarios like UED 12 1. 
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